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ABSTRACT 

In the Local Group (LG), almost all satellite dwarf galaxies that are within the virial radius of the 
Milky Way (MW) and Andromeda (M31) exhibit strong environmental influence. The orbital histories 
of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, 
and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories 
of satellites with M star = 10 3-9 Mg using the ELVIS suite of cosmological zoom-in dissipationless 
simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 
5 — 8 Gyr ago at z = 0.5 — 1. However, they first fell into any host halo typically 7—10 Gyr ago 
at z = 0.7 — 1.5. This difference arises because many satellites experienced “group preprocessing” in 
another host halo, typically of M v [ r ~ 10 10-12 Mg, before falling into the MW/M31 halos. Satellites 
with lower mass and/or those closer to the MW/M31 fell in earlier and are more likely to have 
experienced group preprocessing; half of all satellites with M star < 10 6 M 0 were preprocessed in a 
group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, 
none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during 
reionization (z > 6), and only < 4% were satellites of any other host halo during reionization. Thus, 
effects of cosmic reionization versus host-halo environment on the formation histories of surviving 
dwarf galaxies in the LG occurred at distinct epochs, separated typically by 2 — 4 Gyr, so they are 
separable theoretically and, in principle, observationally. 

Subject headings: cosmology: theory — galaxies: dwarf — galaxies: groups: general — galaxies: 
interactions — Local Group — methods: numerical 


1. INTRODUCTION 

Galaxies in dense environments are more likely to have 
suppressed (quiescent) star-formation rates, more ellipti¬ 
cal/spheroidal morphologies, and less cold gas in/around 
them than galaxies of similar stellar mass, M star , in less 
dense environments. While such environmental effects 
long have been studied for galaxies in massive galaxy 
groups and clusters (for example, Oemler 1974; Dressier 
1980; Dressier & Gunn 1983; Balogh et al. 1997; Blanton 
& Moustakas 2009, for review), the observed effects on 
the dwarf galaxies in the Local Group (LG), in particu¬ 
lar, the satellites within the host halos of the Milky Way 
(MW) and M31, are even stronger (for example, Mateo 
1998; McConnachie 2012; Phillips et al. 2014; Slater & 
Bell 2014; Spekkens et al. 2014). 

Specifically, the galaxies around the Milky Way (MW) 
and Andromeda (M31) show a strikingly sharp transi¬ 
tion in their properties within « 300 kpc, corresponding 
to the virial radii, AW, of the halos of the MW and M31 
for M v i r ss 10 12 Mg (for example, Deason et al. 2012; 
van der Marel et al. 2012; Boylan-Kolchin et al. 2013). 
Within this distance, galaxies transition from (1) having 
irregular to clliptical/spheroidal morphologies, (2) hav¬ 
ing most of their baryonic mass in cold atomic/molecular 
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gas to having little-to-no detectible cold gas, and (3) 
being actively star-forming to quiescent (McConnachie 
2012, and references therein). This environmental tran¬ 
sition of the population is almost complete, with just 
a few exceptions. Four gas-rich, star-forming, irregular 
galaxies persist within the halos of the MW (the LMC 
and SMC) and M31 (LGS 3 and IC 10). However, the 
LMC and SMC are likely on their first infall (Besla et al. 
2007; Kallivayalil et al. 2013), and given their distances 
to M31, LGS 3 and IC 10 may be as well. Furthermore, 4 
- 5 gas-poor, quiescent, spheroidal galaxies exist beyond 
the halos of either the MW and M31: Cetus, Tucana, 
KKR 25, KKs 3 (Karachentsev et al. 2015) and pos¬ 
sibly Andromeda XVIII. While the radial velocities of 
Cetus and Tucana imply that they likely orbited within 
the MW halo (Teyssier et al. 2012), KKR 25 and KKs 
3 are much more distant at ~ 2 Mpc. The fact that 
almost all of the satellite galaxies within the MW/M31 
halos show such strong environmental effects is particu¬ 
larly striking given that, other than KKR 25 and KKs 3, 
all known galaxies at M star < 10 9 Mg that are isolated 
(not within 1500 kpc of a more massive galaxy, and thus 
not strongly influenced by environmental effects) are ac¬ 
tively star-forming (Geha et al. 2012). Thus, the MW 
and M31 halos exert the strongest environmental influ¬ 
ence on their galaxy populations of any observed systems, 
making the LG one of the most compelling laboratories 
to study environmental effects on galaxy evolution. 

Several environmental processes within a host halo reg¬ 
ulate the gas content, star formation, morphology, and 
eventual tidal disruption of satellite galaxies. Gravi¬ 
tationally, the strong tidal forces of the host halo will 
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strip mass from the satellite (subhalo) from the outside- 
in (Dekel et al. 2003; Diemand et al. 2007; Wetzel & 
White 2010). In addition, the dense collection of satel¬ 
lites within a host halo can drive impulsive gravitational 
interactions with each other (Farouki & Shapiro 1981; 
Moore et al. 1998), and satellites can merge with one 
another (Angulo et al. 2009; Wetzel et al. 2009b, a; Dea¬ 
son et al. 2014a). Moreover, tidal shocking and reso¬ 
nant interactions with the host’s galactic disk can lead 
to particularly efficient morphological evolution, coring, 
stripping, and disruption (Mayer et al. 2001; D’Onghia 
et al. 2010; Zolotov et al. 2012). Hydrodynamically, if 
the host halo contains thermalized hot gas, this can strip 
and heat the extended gas from the orbiting satellite sub¬ 
halo (Balogh et al. 2000; McCarthy et al. 2008), leading 
to reduced gas cooling/accretion into the satellite’s disk 
(Larson et al. 1980). More drastically, given a sufficiently 
high density of hot gas and high orbital velocity, ram- 
pressure can strip cold gas directly from the satellite’s 
disk (Gunn & Gott 1972; Abadi et al. 1999; Mayer et al. 
2006; Chung et al. 2009; Tonnesen & Bryan 2009). Fur¬ 
thermore, feedback from stars and/or AGN within the 
satellites can drive galactic winds that can enable these 
environmental process to operate even more efficiently 
(for example, Bahe & McCarthy 2015). 

Understanding the relative efficiency of the above envi¬ 
ronmental processes, including the timescales over which 
they have operated, requires understanding in detail the 
orbital and virial-infall histories of the current satellite 
population in the context of the hierarchical structure 
formation of ACDM. While some authors examined the 
virial-infall times of LG-like satellites in cosmological set¬ 
tings (for example, Lux et al. 2010; Rocha et al. 2012), 
such works used cosmological zoom-in simulations of one 
or two MW-like halos, which does not model the environ¬ 
ment of the MW/M31 pair in the LG or allow for good 
statistics. In addition, hierarchical growth means that 
many satellites may have been in a group before they fell 
into the MW/M31 halos. Environmental processing in 
such groups could help to explain the high efficiency and 
near completeness of environmental effects on the satel¬ 
lite population. Several authors explored the importance 
of this “group preprocessing” on satellites within massive 
groups and clusters (for example, Zabludoff & Mulchaey 
1998; McGee et al. 2009; Hou et al. 2014). In particu¬ 
lar, Wetzel et al. (2013) found that group preprocessing 
alone largely can account for the fact that satellites in 
more massive groups/clusters are more likely to be qui¬ 
escent. However, on mass scales of MW/M31 halos, the 
impact of group preprocessing on dwarf galaxies remains 
largely unexplored. Using a cosmological zoom-in simu¬ 
lation of a single MW-like halo Li & Hclmi (2008) found 
that ~ 1/3 of satellites fell in as part of a group, and us¬ 
ing the two Via Lactea simulations, Slater & Bell (2013) 
similarly found that many satellites are organized into 
small groups with correlated infall. 

If some of the satellites in the MW/M31 halos fell in 
as part of a group, this would have several implications 
for their subsequent evolution and spatial distribution. 
For instance, group infall could have caused many of the 
strong associated in phase-space coordinates between the 
observed satellites (and streams) in the LG (Lynden-Bell 
& Lynden-Bell 1995; D’Onghia & Lake 2008; Klimen- 
towski et al. 2010), including the disk-like configurations 


of satellites around the MW and M31 (for example, Libe- 
skind et al. 2005; Lovell et al. 2011; Fattahi et al. 2013; 
Ibata et al. 2013; Pawlowski & Kroupa 2014). Further¬ 
more, group infall could have driven mergers between 
satellites after infall (for example, Deason et al. 2014a). 

In addition to the above environmental processes that 
operate within a host halo, cosmic reionization may have 
had a lasting impact on formation histories of dwarf 
galaxies in the LG, by heating/removing the gas from 
low-mass halos whose virial temperatures were below 
that of the ultra-violet photoionization background, thus 
quenching star formation in the lowest-mass galaxies and 
leaving them as relics of reionization (Bullock et al. 2000; 
Gnedin 2000). Many ongoing observational efforts aim 
to use the current stellar populations of faint and ultra¬ 
faint satellites in the LG to test the impact that reion¬ 
ization may have had on their star-format ion histories 
at z > 6 (for example, Brown et al. 2014; Weisz et al. 
2014). However, a long-standing challenge for such stud¬ 
ies is whether one can separate the effects of cosmic reion¬ 
ization from those of the host-halo environment on the 
formation histories of surviving satellite galaxies. 

In this work, we examine the orbital and virial-infall 
histories of the current satellite galaxies of the LG in a 
fully cosmological and hierarchical context, including the 
impact of group preprocessing and implications for using 
such galaxies as probes of reionization. Specifically, we 
will address the following questions for the satellites in 
the halos of the MW and M31: 

(a) When did they fall into the MW/M31 halo, and 
when did they first fall into any host halo? 

(b) What fraction were within their MW/M31 halo, 
or any other host halo, during cosmic reionization 
(z > 6)? 

(c) What fraction were in a group prior to falling into 
the MW/M31 halo? 

(d) What role does group infall play in driving mergers 
between satellites? 

2. NUMERICAL METHODS 
2.1. Simulations 

To study the orbital histories of satellite dwarf galaxies, 
we use ELVIS (Exploring the Local Volume in Simula¬ 
tions), a suite of cosmological zoom-in iV-body simula¬ 
tions that are targeted to modeling the LG (Garrison- 
Kimmel et al. 2014). ELVIS was run using GADGET- 
3 and GADGET-2 (Springel 2005), with initial con¬ 
ditions generated using MUSIC (Hahn & Abel 2011), 
all with ACDM cosmology based on WMAP7 (Larson 
et al. 2011): cr 8 = 0.801, Q matter = 0.266, U A = 0.734, 
n s = 0.963 and h = 0.71. 

ELVIS contains 48 dark-matter halos of masses similar 
to the MW or M31 (M v ; r = 1.0 —2.8 x 10 12 M 0 ) within a 
zoom-in volume of radius > 4 J? v ; r of each halo at z = 0. 
Half of these halos are located in zoom-in regions that 
were selected to contain a pair of halos that resemble 
the masses, distance, and relative velocity of the MW- 
M31 pair, while the other half are single isolated halos 
matched in masses to the paired ones. These zoom-in 
regions are selected from a suite of simulations, each a 
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cube with side length 70.4 Mpc. Within the zoom-in re¬ 
gions, the particle mass is 1.9 x 1O 5 M 0 and the Plummer- 
equivalent force softening is 140 pc (comoving at z > 9, 
physical at z < 9). Additionally, three of the isolated ha¬ 
los were re-run at higher resolution, with particle mass 
of 2.4 x 1O 4 M 0 and force softening of 70 pc. Using these 
simulations, we checked that resolution does not signif¬ 
icantly affect any results in this work. See Garrison- 
Kimmel et al. (2014) for more details on ELVIS. 

Throughout this work, unless otherwise stated, we use 
the paired halos, which more accurately capture the en¬ 
vironment, assembly history, and massive satellite pop¬ 
ulation (LMC/M33-like satellites) of the LG. More pre¬ 
cisely, we use 10 halo pairs (20 halos total), ignoring 
the pairs Siegfried & Roy and Serena & Venus because 
they contain a massive halo within 1.2 Mpc that is not 
representative of the LG (Garrison-Kimmel et al. 2014). 
Henceforth, we refer to these paired halos as “MW/M31 
halos”, and we do not further differentiate between the 
MW and M31 given their similar masses (van der Marel 
et al. 2012; Boylan-Kolchin et al. 2013). In the Appendix, 
we compare our results for paired versus isolated halos. 

2.2. Finding and tracking (sub)halos 

ELVIS identifies dark-matter (sub)halos using the 
six-dimensional halo finder ROCKSTAR (Behroozi 
et al. 2013b) and constructs merger trees using the 
CONSISTENT-TREES algorithm (Behroozi et al. 2013c). 
For each halo that is not a subhalo (see below), we assign 
a virial mass, M v ; r , and radius, R V1I , using the evolution 
of the virial relation from Bryan & Norman (1998) for 
our ACDM cosmology. At z = 0, this corresponds to an 
overdensity of A critical = 97 (A matt er = 363) times the 
critical (matter) density, while at z > 3 it asymptotes to 
A C ritical ~ A ma tter ~ 178. 

We define a “host halo” as an isolated halo that can 
host (lower-mass) subhalos within in, and a “subhalo” 
as a halo whose center is inside R v i r of a (more massive) 
host halo. When a (sub)halo passes within R v i r of a host 
halo, the (sub)halo becomes its “satellite” and experi¬ 
ences “virial infall”. 

For each (sub)halo, we assign its primary progenitor 
(main branch) as the progenitor that contains the most 
total mass summed from the (sub)halo masses over all 
preceding snapshots in that branch. We then compute 
the peak mass, M pea jj, as the maximum instantaneous 
mass that a (sub)halo ever reaches along the history of 
its primary progenitor. For subhalos, M pea k almost al¬ 
ways occurs before virial infall. As explored in Garrison- 
Kimmel et al. (2014), the numerical resolution of ELVIS 
does not significantly affect its (sub)halo catalogs at (and 
even below) M pea k > 10 s M 0 , the limit in this work. 

2.3. Assigning stellar mass to subhalos 

Our goal is to map luminous galaxies to the dark- 
matter subhalos in ELVIS. The relation between stel¬ 
lar mass and subhalo mass (or maximum circular veloc¬ 
ity) for dwarf galaxies is highly uncertain, likely with 
significant scatter, especially for our lowest-mass subha¬ 
los, some of which might not host any luminous galaxies, 
a manifestation of the long-standing “missing satellites 
problem” (Klypin et al. 1999). Nonetheless, we use the 
relation from abundance matching to ELVIS subhalos in 


Garrison-Kimmel et al. (2014), which is based on that of 
Behroozi et al. (2013a) but is modified at the low-mass 
end according to the observed stellar-mass function of 
Baldry et al. (2012). At the mass scales of dwarf galax¬ 
ies, this leads to M star oc This modification re¬ 

produces the observed mass function at M star < 10 9 M 0 
in the LG, especially if one accounts for observational 
incompleteness (Tollerud et al. 2008; Hargis et al. 2014). 
However, we present most results as a function of both 
M star and Af pea k, given the uncertainties of abundance 
matching at these low masses. For reference, each of 
the MW/M31 halos hosts an average of 230, 28, and 3 
satellites with M star = 10 3-5 , 10 5-7 , and 10 7-9 M 0 , re¬ 
spectively, at z = 0 (this is similar for the paired and 
isolated halos). 

In selecting the stellar-mass ratio for defining “major” 
groups or mergers in the histories of satellites in Sec¬ 
tions 4.1 and 5, we assume that the slope (but not nec¬ 
essarily normalization) of this relation does not evolve, 
motivated by the lack of strong evolution observed for 
slightly more massive galaxies (for example, Leauthaud 
et al. 2012; Hudson et al. 2013) and the lack of obser¬ 
vational evidence to suggest otherwise. We define major 
mergers as those for which the ratio in M star >0.1. This 
broadly corresponds to mass ratios at which the lower- 
mass companion likely has significant dynamical effect 
on the more massive galaxy (for example, Hopkins et al. 
2010; Helmi et al. 2012; Yozin & Bekki 2012) and for 
which recent mergers are likely to be observable. Given 
our relation between M star and M pea k, this corresponds 
to a ratio in M pea k > 0.3. 

3. VIRIAL-INFALL TIMES OF SATELLITES 

We start by investigating the virial-infall times of satel¬ 
lite dwarf galaxies at z = 0, to understand how long they 
have been satellites within a host halo. This timescale 
has several important implications. First, it provides in¬ 
sight into how long satellites have experienced environ¬ 
mental processes that cause the observed depletion of 
gas, quenching of star formation, transition of morphol¬ 
ogy, and (potentially) stripping of stars. Second, it tells 
us what fraction of satellites at z = 0 were satellites 
within a host halo during the epoch of cosmic reioniza¬ 
tion. This allow us to understand whether the differ¬ 
ing effects of reionization versus host-halo environment 
on surviving satellites occurred at distinct epochs in the 
formation histories of surviving satellites. 

3.1. Determining virial infall for satellites 

While environmental processes clearly affect satellite 
galaxies within R v i r of the MW/M31 halos, whether such 
environmental processing occurs in lower-mass host ha¬ 
los (M v i r <C 10 12 M 0 ) remains unclear. Thus, we in¬ 
vestigate two metrics of virial infall. First, we examine 
“first infall”: when a satellite first became a satellite 
within any host halo, that is, first crossed within the 
virial radius of any halo more massive than itself. We re¬ 
fer to this redshift as Zfirstinfa.il or time as tfi rs t infall, with 
^firstlnf a ii = ^now - tfirst mfaii • We also examine “MW/M31 
infall”: when a satellite first became a satellite in its 
host MW/M31 halo. We refer to this as z MW /M3i infall or 

^MW/M31 infall; with f]y[\Y/M31 infall = ^now f MW/M3I infall■ 
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3.2. Dependence of virial-infall time on satellite mass 

We first examine how the virial-infall times of satel¬ 
lites at z = 0 depend on their mass. Figure 1 shows 

i&stlnfall ( left ) and 4 n W/M31mfall ( ri g ht )> OT ^first infall and 

^MW/M 3 i infail on the right axes, as a function of satellite 
Mstar, or subhalo M pea k on the top axes. 

For both virial-infall metrics, lower-mass satellites fell 
in earlier, though significant scatter persists at all masses. 
This trend with M pea k (or M sta r) is a natural result of 
hierarchical structure formation, for two reasons. First, 
halos of a given Af pea k are more common at later cosmic 
time. Thus, higher-mass satellites are more likely to have 
formed, and subsequently fallen in, at later time. Sec¬ 
ond, satellites with higher M pea k have shorter dynamical- 
friction lifetimes (for fixed host-halo mass) before they 
tidally disrupt or merge with the host (Boylan-Kolchin 
et al. 2008; Jiang et al. 2008; Wetzel & White 2010). 

Our lowest-mass (ultra-faint) satellites first fell into 
any host halo typically ~ 10 Gyr ago at 2 ~ 1.7, and they 
first fell into the MW/M31 halo ~ 7.5 Gyr ago at z ~ i. 
By contrast, our highest-mass satellites (corresponding 
to the LMC, SMC, NGC 205, M32) have tf( r n s c t e infall ~ 
6.5 Gyr (z first in f a ii ~ 0.8) and M 3 iinf a ii ~ 5 G y r 
(^mw/M 3 i inf a ii ~ 0.5). Thus, satellites first fell into any 
host halo a few Gyr before they fell into the MW/M31 
halo, as we will examine further in Section 4.3. This 
result is generic for hierarchical structure formation. 

Overall, satellite dwarf galaxies at z = 0 typically have 
evolved as a satellite in a host halo for over half of their 
entire history, so the host-halo environment typically has 
had significant time to affect their evolution. 

3.3. Dependence of virial-infall time on satellite 
distance 

We next explore how the above virial-infall times of 
satellites depend on their current distance from the 
MW/M31 host. Previous analyses at higher mass scales 
in cosmological simulations showed that satellites at 
smaller distances tend to have fallen in earlier (for exam¬ 
ple, Gao et al. 2004; Smith et al. 2012; Oman et al. 2013). 
This is because (1) R v i r of the host halo was smaller 
at earlier time, and (2) satellite orbits get dragged to 
smaller distance over time via dynamical friction. Such 
a trend for satellites in the LG would be important for 
several reasons. First, ultra-faint (M star < 10 5 Mg) 
satellites currently are observable only within the inner 
~ 50 kpc of the MW halo, so it is possible that they 
had systematically earlier infall times than the median 
in Figure 1. Second, such a correlation with distance 
would provide a statistical proxy for an environmental 
evolutionary sequence (since the time of infall) for the 
observable satellite population. 

Figure 2 shows ^" c t e infaU (top row) and 4^ e /M31infaU 
(bottom row) as a function of distance to the host’s cen¬ 
ter, d , as scaled to the host’s virial radius, R vlr , at z = 0. 
We compute these quantities in bins of d/R vn - for each 
MW/M31 halo, and using that the median R vlr in ELVIS 
is 300 kpc, we also show the (approximate) dependence 
on d along the top axes. For both infall metrics and 
at all masses, satellites closer to the host fell in ear¬ 
lier, though with significant scatter. Overall, the trend 
for MW/M31 infall is slightly stronger than for first in¬ 


fall, as expected given that we measure d with respect 
to the MW/M31 center. The gradient in and 

*MW/M 3 i infaii fr° m d/R vu - = 0 to 1 for massive satellites 
is 6.5 and 4.5 Gyr, respectively, while for our lowest- 
mass satellites it is 2.5 and 4 Gyr. Thus, the corre¬ 
lation of infall time with distance is stronger for more 
massive satellites, because (1) they fell in more recently, 
making them less smeared out in orbital phase space, 
and (2) they experience more efficient dynamical fric¬ 
tion. For our lowest-mass (ultra-faint) satellites, those 
in our smallest distance bin, where they are observable, 
experienced first infall typically ~ 11 Gyr ago at z ~ 2.2, 
and they first fell into the MW/M31 halo ~ 9 Gyr ago 
at z ~ 1.5, so these are slightly earlier than Figure 1. 

We conclude that a satellite’s distance does provide a 
statistical proxy for its virial-infall time, and therefore, 
the distribution of distances for the satellite population 
at z = 0 provides a proxy for an environmental evolu¬ 
tionary sequence, especially for more massive satellites. 
However, note one caveat: despite the correlation with 
distance, satellites currently near 7? v ; r have been satel¬ 
lites for quite a while, typically 3 — 8 Gyr, depending on 
mass. Thus, the satellite population near l? v i r is not only 
a recent-infall population, but rather, it is a superposi¬ 
tion of inward- and outward-orbiting satellites, some of 
which already experienced a pericentric passage. Note 
that any given satellite spends most of its orbital time 
near apocenter. (For reference, the virial crossing time, 
Rvir/Vvii, is « 2 Gyr at z = 0.) 

3.4. Implications for dwarf galaxies during cosmic 
reionization 

The virial-infall times in Figures 1 and 2 have impor¬ 
tant implications for understanding the relative effects of 
cosmic reionization versus host-halo environment on sur¬ 
viving dwarf galaxies, in particular, whether the effects 
of these two processes occurred at distinct epochs during 
the formation histories of these galaxies. In these figures, 
the dashed line at z = 6 represents the end of cosmic 
reionization as constrained by various observations (for 
example, Robertson et al. 2013, and references therein). 
Across all masses and distances, none of the satellites at 
z = 0 were within R v j r of their MW/M31 halo any time 
during reionization. Furthermore, < 4% were within i? v j r 
of any host halo during reionization. There are regimes 
where this fraction is somewhat higher, the highest being 
10% for satellites at M star = 10 5-8 M 0 and d/R v i r < 0.1, 
but this fraction is typically only a few percent across our 
range of mass and distance. Thus, essentially none of the 
dwarf galaxies in the LG were within f? v i r of a host halo 
during reionization, such that they would have experi¬ 
enced strong environmental effects at that time. 

To understand this result in more detail, we also exam¬ 
ine how close the dwarf galaxies came to a more massive 
halo during reionization. Thus, we select all satellites in 
the MW/M31 halos at z = 0 and trace them back to 
z > 6, when almost all such dwarfs were isolated (non¬ 
satellite) halos. (We are able to track all (sub)halos back 
to z > 6, except for 6% of those at M star < 10 4 M 0 , 
which formed after z = 6.) At all z > 6 (for which ELVIS 
contains 6 snapshots), we then compute the nearest dis¬ 
tance, d n earest, that each satellite’s progenitor came to 
the center of any neighboring halo that is more massive 
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Figure 1. Time since virial infall for satellite dwarf galaxies at z = 0 as a function of their stellar mass, M s t ar , or subhalo peak mass, 
M pea k: time since first crossing within R v ; r of any host halo, thus including group preprocessing (left), or time since first crossing within 
7i vir of the MW/M31 halo (right). Curves show median, shaded regions show 68%, 95%, and 99.7% of the distribution. Satellites at z = 0 
have been satellites typically for over half of their history. Lower-mass satellites fell in earlier, though with large scatter. Dotted line at 
2 = 6 marks the end of cosmic reionization; during reionization, < 4% of current satellite galaxies were a satellite in a host halo, and none 
were in the MW/M31 halo, demonstrating the effects of reionization versus host-halo environment occurred at distinct epochs, separated 
typically by 2 — 4 Gyr, and thus are separable in time theoretically and, in principle, observationally, during satellites’ evolutionary histories. 
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Figure 2. Time since virial infall for satellite dwarf galaxies at z = 0 as a function of their distance to their host, d, or as scaled to the 
host’s virial radius, R V1T , at z = 0: time since first crossing within R v [ r (z ) of any host halo (top row) or within R v [ r (z) of the MW/M31 
halo (bottom row). Left-to-right columns show satellites of increasing M s tar- Curves shows the median, shaded regions show 68%, 95%, 
and 99.7% of the distribution. At all masses, satellites closer to the host fell in earlier, though with significant scatter. Dotted line at 
z = 6 indicates the end of cosmic reionization; even for the lowest-mass (ultra-faint) satellites at the smallest distances, where they are 
observable, < 4% were a satellite in any host halo during reionization, and none were in the MW/M31 halo during reionization. 
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and thus feasibly could induce environmental effects. 

Figure 3 shows the cumulative distribution of ^nearest 
in comoving units (left) and scaled to R v i r of the near¬ 
est, more massive halo (right). The thick solid curve 
shows the average across the paired MW/M31 halos. We 
show the average for all satellites across our mass range, 
Af star = 10 3-9 M 0 , as we find no significant dependence 
on satellite mass. The typical ^nearest at z > 6 was large: 
3 Mpc comoving (~ 400 kpc physical), or 500 R v i r . More¬ 
over, only ss 5% of these galaxies came within 1 Mpc 
comoving, or 100f? v ; r , of a more massive halo. The thin 
solid curves show each MW/M31 pair, highlighting the 
factor of ~ 2 scatter in these distributions. For reference, 
the thick dotted curve shows the average for the isolated 
MW/M31 halos, which we discuss in the Appendix. 

If strong environmental effects (such as star-formation 
quenching) on dwarf galaxies are confined to within R v - lr 
(or even larger) of a host halo at these redshifts, then 
the results of Figures 1,2, and 3 indicate that such en¬ 
vironmental processing occurred only at z < 6 during 
the histories of satellites at z = 0, and more typically, 
at z < 3 (for ss 84% of all satellites). The properties of 
the LG support this result, because a strong transition 
in morphology, star formation, and gas content of dwarf 
galaxies, as induced by the host-halo, occurs only within 
~ 300 kpc (~ R v h ) of the MW or M31. 

Given that cosmic reionization ended by z = 6, we thus 
conclude that the effects of reionization occurred well 
before those of the host-halo environment for surviving 
satellites. Specifically, for surviving satellites, the dura¬ 
tion between the end of reionization and first crossing 
within R vil of a host halo was typically > 2 Gyr for first 
infall into any host halo (including group preprocessing) 
and > 4 Gyr for first infall into the MW/M31 halo. This 
result strongly motivates the use of faint and ultra-faint 
satellites as probes of reionization, for example, by mea¬ 
suring the (potential) effect of reionization on their star- 
formation histories as derived from their current stellar 
populations (for example, Brown et al. 2014; Weisz et al. 
2014), if such star-formation histories can demonstrate 
that any individual satellite quenched at z > 6 or that 
any common features (such as quenching) occurred at 
£ > 3 for a population of several satellites (given that 
~ 84% surviving satellites fell into a host halo at z < 3). 

4. GROUP INFALL AND PREPROCESSING 

In the previous section, we showed that many satel¬ 
lite galaxies first became satellites significantly prior to 
falling into the MW/M31 halo. Thus, many satellites 
spent significant time in another host halo, which may 
have environmentally “preprocessed” them prior to their 
joining the MW/M31 halo. We now explore what frac¬ 
tion of all satellites at z = 0 were preprocessed as a 
satellite in another host halo (group) prior to falling into 
the MW/M31 halo. We examine two such metrics: the 
fraction of all current satellites that were a satellite in 
another host halo (1) any time before falling into the 
MW/M31 halo, or (2) at the time of falling into the 
MW/M31 halo. The difference between these is driven 
by “ejected” or “backsplasli” satellites that fell into an¬ 
other host halo and then orbited out beyond its f? v ir 
before falling inside R V1I of the MW/M31 halo. Recent 
work indicates that these satellites are affected environ¬ 
mentally in similar ways as satellites that remain within 


R v ir (for example, Ludlow et al. 2009; Knebe et al. 2011; 
Teyssier et al. 2012; Bahe et al. 2013; Wetzel et al. 2014). 
Thus, we examine both preprocessing metrics. 

4.1. Dependence of group preprocessing on satellite 

mass 

Figure 4 (left) shows both of the above preprocessed 
fractions as a function of satellite M sta r or M pea k. First, 
the red short-dashed curve shows the fraction of satellites 
that were a satellite in a group at the time of falling into 
the MW/M31 halo. For our highest-mass satellites, this 
fraction is relatively low (~ 10%) but increases signifi¬ 
cantly to 30% for our lowest-mass satellites. Thus, ~ 1/3 
of all faint and ultra-faint satellites were in a group when 
they fell into the MW/M31 halo, in good agreement with 
the results of Li & Helmi (2008), which was based on a 
cosmological simulation of a single MW-like halo. 

Second, the blue long-dashed curve shows the fraction 
of satellites that were a satellite in a group any time 
before falling into the MW/M31 halo. This fraction is 
significantly (~ 2x) higher, because of the large fraction 
of satellites whose orbits brought them beyond R vir of 
their preprocessing group. For our higliest-mass satel¬ 
lites, ~ 30% were preprocessed by a group before joining 
the MW/M31 halo. Again, this preprocessed fraction 
increases at lower mass, being ~ 60% for faint and ultra¬ 
faint satellites. Thus, half of all satellite galaxies with 
Af star < 10 6 M 0 were preprocessed as a satellite in a 
group prior to joining the MW/M31 halo. 

We also explore how many current satellites in the 
MW/M31 halo were the central (most massive) galaxy 
in such an infalling group. That is, we identify satellites 
in the MW/M31 halo at z = 0 that hosted their own 
major satellite(s) when they fell into the MW/M31 halo. 
(By “major”, we mean that M star differs by less than a 
factor of 10, as detailed in Section 2.3). Figure 4 (right) 
shows this fraction via the green long-short-dashed curve, 
which is 5 — 10% across our mass range. This fraction 
increases only weakly with mass because of the combi¬ 
nation of (1) a nearly mass-independent distribution of 
Alpe a k, s a teiiite/Mpe a k, host bi ACDM, (2) our power-law 
Af star — Afp ea k relation, and (3) our requirement that 

A/stax, satellite/Af s t a r, host A 0.1. 

For comparison, Figure 4 (right) also shows the same 
red short-dashed curve as in the left panel and the black 
dot-dashed curve shows the sum of the two curves, in¬ 
dicating the total fraction of satellites that fell into the 
MW/M31 halo as part of a major group. This overall 
fraction is significant at 20 — 40% across our mass range. 

4.2. Dependence of group preprocessing on satellite 

distance 

We next explore how the above preprocessed fractions 
vary with the current distance of satellites from their 
host. The two panels in Figure 5 shows the same two pre- 
processed fractions as in Figure 4 (left), but as a function 
of d/Rv ir at 2 = 0, similar to Figure 2. Again, we com¬ 
pute these quantities in bins of d/R v i r for each MW/M31 
halo, though we also show the dependence on d (using 
the median R v ; r across the MW/M31 halos) along the 
top axes. At any distance, lower-mass satellites are more 
likely to have been preprocessed. Moreover, at nearly 
all satellite masses, those closer to the host center are 
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Figure 3. For all satellite galaxies at z = 0 with M s tar = 10 3-9 Mq, cumulative distribution of the distance to the nearest, more massive 
halo, ^nearest, that they experienced during cosmic reionization (z > 6). Left panel shows comoving distance, and right panel shows this 
distance scaled to R v i r of the nearest halo. Solid thick curve shows average over all satellites in the paired MW/M31 halos, while thin 
curves shows satellites in each pair, to indicate the pair-to-pair scatter. We find no dependence on satellite mass. The typical distance was 
3 Mpc comoving (~ 400 kpc physical), or 500-R v ir5 at these distances, dwarf galaxies in the Local Group do not show strong environmental 
influence. These results strongly support that the effects of the host-halo environment occurred well after the effects of reionization during 
the histories of surviving satellites. For comparison, dotted curve shows the average across the isolated MW/M31 halos, whose satellites 
experienced ~ 1/2 the distance (see Appendix). 
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Figure 4. Fraction of all satellite dwarf galaxies at z = 0 that experienced various aspects of group preprocessing prior to falling into the 
current MW/M31 halo, as a function of satellite stellar mass, M s tar, or subhalo peak mass, M pea k. Left: Fraction that were a satellite in 
a group any time prior to (blue long-dashed) or at the time of (red short-dashed) falling into the MW/M31 halo. Curves show average 
over the paired MW/M31 halos, while shaded regions show standard deviation from halo-to-halo scatter. The difference between the two 
curves is driven by ejected/backsplash satellites that were once within another host halo but then orbited beyond its R v i r . Lower-mass 
satellites are more likely to have been preprocessed in a group: over half of the lowest-mass (faint and ultra-faint) satellites were. Right: 
Red short-dashed curve shows same as in the left panel, while green long-short-dashed curve shows the fraction that were the central (most 
massive) galaxy in a group that contained a major satellite (at least O.lx its M s t ar ) at the time of falling into the MW/M31 halo. Black 
dot-dashed curve show the sum of the two curves. 
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Figure 5. Fraction of all satellite dwarf galaxies at z = 0 that were a satellite in another host halo any time prior to (left) or at the time 
of (right) falling into the MW/M31 halo, as a function of their current distance to the host, d, or as scaled to the host’s virial radius, i? v ir- 
Curves show average over the paired MW/M31 halos in bins of satellite M s tar, and error bars show 68% uncertainty in this fraction for a 
beta distribution (the halo-to-halo scatter is similar to Figure 4). For both group-preprocessing metrics, and at across all masses, satelltes 
closer to the host are more likely to have been preprocessed. 
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more likely to have been preprocessed, with a nearly 2x 
increase from d/R v i r = 1 to 0.1 for our lowest-mass satel¬ 
lites. Most likely, this gradient arises because a satellite 
that fell in as part of a group remained bound to the 
(more massive) group for some time after infall, so it ex¬ 
perienced more efficient dynamical friction that dragged 
it to the center of the MW/M31 halo more rapidly. 

Overall, given that the faintest galaxies are observable 
only at small distances within the MW halo, most likely 
about half of them were a satellite in another group (s) 
before / during infall into the MW halo. 

4.3. Duration and host-halo mass of group preprocessing 

We next examine in more detail the durations and 
host-halo masses that satellites at z = 0 experienced 
during their group preprocessing, in order to understand 
better its potential impact on their evolution. 

Figure 6 (top) shows the distribution of maximum 
host-halo mass that satellites experienced during group 
preprocessing as a function of their M star - The typi¬ 
cal preprocessing group had M v j r ~ 10 11 M 0 , with 68% 
spread of lO 10- ^ 2 M 0 , largely independent of satellite 
mass, though with scatter to M v - lr at lower M star - Com¬ 
bined with Figure 4 (left), this means that ~ 25% of 
all satellites at z = 0 were preprocessed in groups of 
M v i r > 10 11 M 0 , comparable to the LMC or M33. 

The total population of preprocessed satellites at z = 0 
originated from typically 30 infalling groups, with most 
groups containing 2 — 5 satellites. However, about half 
of all preprocessed satellites fell in via the 2 — 3 most 
massive groups (M v i r > 10 11 M 0 , similar to the LMC 
or M33), and the most massive infalling group typically 
brought in ~ 25 — 30% of all preprocessed satellites, or 
~ 13% of all satellites in total. 

However, for half of the preprocessed satellites, their 
preprocessing host does not survive to z = 0 but in¬ 
stead merges/disrupts within the MW/M31 halo. Thus, 
surviving preprocessed satellites do not necessarily have 


their preprocessing host galaxy surviving near them at 
z = 0; instead, the host galaxy could persist in a dis¬ 
rupted stream configuration. On the other hand, par¬ 
ticularly massive (M pea k > 10 11 M 0 , comparable to 
the LMC) satellites that do survive to z = 0 typically 
brought in ~ 7% of overall satellite dwarf population in 
the MW/M31 halo (Deason et al. 2015). 

Satellites can be preprocessed by the other MW/M31 
halo in the pair, if they fell into one of the MW/M31 ha¬ 
los, orbited beyond its R v i r , and then fell into the other 
MW/M31 halo. However, this accounts for < 2% of all 
preprocessed satellites across our mass range, in general 
agreement with Knebe et al. (2011), so this is not a par¬ 
ticularly important population. 

Figure 6 (bottom) shows the distribution of the du¬ 
ration that satellites spent in a preprocessing host halo. 
The typical duration was ~ 1.2 Gyr, with 68% spread of 
0.5-3.5 Gyr. At M star > 10' M 0 , no satellites were pre- 
processed longer than 1.8 Gyr, while at M star < 10 7 M 0 , 
the scatter increases significantly, with some satellites 
having experienced up to 7 Gyr of preprocessing. The 
small durations at M stal > 10 7 M 0 likely arise because 
those satellites have M pea k that approaches that of their 
preprocessing host halo, corresponding to shorter dy¬ 
namical friction lifetimes, so they could not have been 
preprocessed too long without merging/disrupting. 

Overall, most preprocessing occurred within groups of 
M v j r = 10 10-12 M 0 , masses that feasibly could influence 
satellite galaxies, though environmental effects at these 
masses remain poorly understood. Furthermore, the typ¬ 
ical preprocessing duration was 0.5 — 3.5 Gyr, compara¬ 
ble to typical timescales over which such satellite dwarf 
galaxies are quenched environmentally (Fillingham et al. 
2015; Wetzel et al. 2015). Thus, we conclude that such 
group preprocessing before joining the MW/M31 halo is 
likely an important component in the evolution of satel¬ 
lite dwarf galaxies. 
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Figure 6. For satellite dwarf galaxies that were preprocessed in 
a group before falling into the MW/M31 halo, the maximum host- 
halo mass experienced during preprocessing (top) and the total 
duration of preprocessing (bottom), as a function of satellite stellar 
mass, Mstar, or subhalo peak mass, M pea k. Solid curves shows 
median, shaded regions show 68, 95, 99.7% of the distribution. 
Most preprocessing occurred in groups with M v j r = 10 10—12 Mq, 
masses at which environmental effects are feasible, though poorly 
understood. The typical group preprocessing duration was 0.5 — 
3.5 Gyr, though some satellites at M s tar < 10 7 Mq experienced 
significantly longer duration. 

5. GROUP INFALL DRIVES SATELLITE-SATELLITE 
MERGERS 
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Figure 7. For all satellite galaxies at z = 0 with Mstar = 
1Q3—9 Mq that experienced a major merger with another satellite 
after falling into the MW/M31 halo, the fraction of such mergers 
that fell in as part of the same host halo (group), as a function 
of the time since (or redshift of) the last major merger, Tlmm 
( or zlmm)- We do not find any significant dependence on satel¬ 
lite mass. Infalling groups drive the majority (60 — 90%) of all 
satellite-satellite mergers. 



satellite-satellite mergers. 

Figure 7 shows, for all satellites that experienced a 
major merger after falling into the MW/M31 halo, what 
fraction occurred between galaxies that fell into the 
MW/M31 halo in the same group, as a function of the 
time since the last major merger, Tlmm- To address 
the limited statistics of major mergers, unlike in the 
rest of this paper, we combine all 48 paired and isolated 
MW/M31 halos in ELVIS, and we bin all satellites across 
Mstar = 10 3 - 9 Mq, given that we do not find any signif¬ 
icant dependence on satellite mass. 

The majority (60 — 90%) of all satellite-satellite ma¬ 
jor mergers occurred between two galaxies that were in 
the same group when they fell into the MW/M31 halo. 
This group infall drives a lower fraction of mergers at 
later cosmic time. This likely relates to the larger de¬ 
lay time between MW/M31 infall and merging at later 
cosmic time (Deason et al. 2014a), which suggests that 
satellites had more opportunity to experience a “chance” 
merger with another satellite in the MW/M31 halo at 
later cosmic time. 

Overall, in addition to driving environmental prepro¬ 
cessing, group infall is also an important catalyst for ma¬ 
jor mergers between satellite galaxies after they fall into 
the MW/M31 halo. 


In Deason et al. (2014a), we showed that most (> 70%) 
satellites in the LG experienced a major merger (M star 
ratio > 0.1) during their history. While most mergers 
occurred prior to falling into the MW/M31 halo, a sig¬ 
nificant fraction were a satellite-satellite merger after in¬ 
fall. Almost all of the latter occurred soon after falling in 
the MW/M31 halo (see Figure 3 in Deason et al. 2014a), 
which suggests that such mergers occurred between satel¬ 
lites with correlated infall histories (for example, Li & 
Helrni 2008; Angulo et al. 2009; Wetzel et al. 2009b, a), 
in particular, that were part of the same group at infall. 
We now demonstrate that group infall drives most such 


6. SUMMARY AND DISCUSSION 
6.1. Summary 

Using the ELVIS suite of cosmological zoom-in dissi¬ 
pationless simulations, we examined the virial-infall his¬ 
tories and group preprocessing of satellites across the 
observed range of masses for dwarf galaxies: M star = 
10 3-9 Mq. While we examined all 48 MW/M31 halos in 
ELVIS, we focused on the 10 halo pairs (20 halos total) 
that most resemble the LG, thus providing good statistics 
in a realistic cosmic setting. We summarize our primary 
results as follows. 
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(a) Virial-infall histories: satellites at z = 0 fell into 
the MW/M31 halo typically 5 — 8 Gyr ago at 
z = 0.5 — 1, though they first fell into any host 
halo typically 7—10 Gyr ago at 2 = 0.7 — 1.5. The 
difference between these infall times arises because 
of group preprocessing. Satellites at lower mass or 
closer to the center of the MW/M31 experienced 
earlier infall times. The latter means that the dis¬ 
tribution of distances of satellites at 2 = 0 provides 
a statistical proxy for an environmental evolution¬ 
ary sequence after infall. Overall, current satellites 
have evolved as satellites within a host halo for over 
half of their entire history. 

(b) Group preprocesssing: a large fraction of satel¬ 
lites were a satellite in a group (another host 
halo), typically of M v ; r ~ 10 10-12 M 0 , for a du¬ 
ration of 0.5 — 3.5 Gyr, before falling into the 
MW/M31 halo. This group preprocessing is es¬ 
pecially common among faint and ultra-faint satel¬ 
lites: at M star < 10 6 M 0 , « 30% of all satellites at 
2 = 0 fell into the MW/M31 halo as a satellite in 
a group, and half of all satellites at 2 = 0 were in 
a group any time before falling into the MW/M31 
halo. Thus, ~ 25% of all satellites at 2 = 0 were 
preprocessed in groups with M v ; r > 10 11 M 0 , com¬ 
parable to the LMC. Satellites closer to the center 
of the MW/M31 are more likely to have experi¬ 
enced group preprocessing. 

(c) Satellite-satellite mergers: group infall drives most 
(60 - 90%) of the satellite-satellite major mergers 
that occurred after falling into the MW/M31 halos, 
as we explored in Deason et al. (2014a). 

(d) Cosmic reionization: none of the surviving satel¬ 
lite dwarf galaxies were within their MW/M31 halo 
during reionization (z > 6), and only < 4% were 
within the virial radius of any host halo during 
reionization. Furthermore, the typical distance to 
the nearest, more massive halo at 2 > 6 was 3 Mpc 
comoving (~ 400 kpc physical), or 500i? v i r . Thus, 
the effects of cosmic reionization versus host-halo 
environment on the formation histories of surviv¬ 
ing dwarf galaxies in the LG occurred at distinct 
epochs separated typically by 2 — 4 Gyr, and are 
separable in time theoretically and, in principle, 
observationally. 

6.2. Discussion 

6.2.1. Impact of group preprocessing on the evolution of 
dwarf galaxies 

The significant fraction of satellite dwarf galaxies that 
experienced group preprocessing may help to explain the 
near-unity fraction of observed satellites in the LG that 
show strong environmental influence: spheroidal mor¬ 
phology, little-to-no cold gas, and quiescent star forma¬ 
tion. However, this depends on the extent to which the 
environmental processes, described in the Introduction, 
operate within low-mass groups of M v j r -C 10 12 M 0 . If 
there is a lower limit in virial mass below which host halos 
do not significantly affect their satellites, then group pre¬ 
processing, even if common, many not be a particularly 
important regulator of the evolution of dwarf galaxies. 


Because such low-mass groups necessarily are composed 
of faint galaxies, few observational works have probed 
the detailed properties of the satellites of such systems, 
though there are now ongoing efforts (for example, Stier- 
walt et al. 2014). Additionally, few theoretical works 
have examined the details of environmental effects in 
such lower-mass host halos, for example, in hydrody¬ 
namic simulations. Based on our results, these would 
be fruitful areas for future investigation. 

6.2.2. Implications for observed associations in the Local 

Group 

The significant group-infall fractions (30 - 60 %) that 
we found may help to explain the many observed asso¬ 
ciations between satellite galaxies (and stellar streams) 
within the halos of the MW and M31. Li & Helrni (2008) 
showed that infalling groups can remain coherent and 
share similar orbital planes for up to ~ 8 Gyr (see also, 
Klimentowski et al. 2010; Sales et al. 2011; Slater & Bell 
2013), which we find is also the typical time that satel¬ 
lites have been within the MW/M31 halo. See also Dea¬ 
son et al. (2015), for detailed phase-space distributions 
at 2 = 0 of infalling LMC-mass groups. 

Observational evidence for associations between dwarf 
galaxies in the LG dates back to Lynden-Bell (1982), 
who demonstrated that the satellites in the MW halo 
appear situated along two great “streams”: the Magel¬ 
lanic stream and the Fornax-Leo-Sculptor stream. The 
discovery of faint and ultra-faint (L < 1O 5 L 0 ) galaxies 
around the MW (Willman et al. 2005; Belokurov et al. 
2006b, 2007a) led to further evidence for such galactic 
associations. For example, the two ultra-faint galaxies 
Leo IV and Leo V are separated by ~ 3 degrees on 
the sky with small offsets in both distance and veloc¬ 
ity (see Fattahi et al. 2013, for more discussion of such 
pairs). Moreover, many works continue to explore (and 
debate) the presence of a planar, disk-like distribution of 
satellites around both the MW and M31 (for example, 
Libeskind et al. 2005; Deason et al. 2011; Lovell et al. 
2011; Fattahi et al. 2013; Ibata et al. 2013; Pawlowski 
& Kroupa 2014), which could be the result of one or 
more infalling groups. Note, however, that Pawlowski & 
McGaugh (2014) examined the satellite configurations in 
ELVIS and found no polar structures analagous to that 
around the MW. Our results also fully support the like¬ 
lihood that the SMC was in a group with the LMC prior 
to MW infall (Kallivayalil et al. 2013). 

In addition to associations between galaxies, many au¬ 
thors have noted observed association between galaxies, 
stellar streams, and/or structures in the stellar halo. For 
example, Newberg et al. (2010) showed that the Orphan 
stellar stream (Grillmair 2006; Belokurov et al. 2007b) 
has a similar distance, velocity, and position as Segue 1. 
Likewise, the proximity of Bootes II in both position and 
velocity to the Sagittarius stream led Koch et al. (2009) 
to suggest that Bootes II may have been stripped from 
the more massive Sagittarius dwarf galaxy. In addition, 
Deason et al. (2014b) (see also Belokurov et al. 2009) 
showed that Segue 2, which is perhaps the least-massive 
known galaxy (Kirby et al. 2013a), is likely associated 
with the large, metal-rich Triangulum-Andromeda over¬ 
density, and Kirby et al. (2013) showed that Segue 2 lies 
off of the tight mass-metallicity relation for most dwarf 
galaxies, which may indicate of group infall. We reiterate 
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that in about half of our cases of group infall, the lower- 
mass satellite from the group survives to z = 0 in the 
MW/M31 halo, while the more massive (primary) galaxy 
merges/disrupts, which could lead to such observable as¬ 
sociations between (surviving) galaxies and (disrupted) 
stellar streams. Overall, the enhanced evidence for as¬ 
sociations between the lowest-mass satellites agrees well 
with our predictions (Figure 4). 

6.2.3. Effects of numerical resolution and baryonic physics 

For any simulation, the survivable lifetime of a satellite 
depends on how well the simulation resolves it, which can 
lead to prematurely merging/disrupting (as compared 
with a real galaxy), especially for lower-mass satellites 
and/or those that fell in at high redshift. If this were 
a strong numerical effect in ELVIS, then the virial-infall 
times for surviving satellites at z = 0 would be biased 
to lower redshifts. However, using the three isolated 
MW/M31 halos in ELVIS that were re-run at 2x higher 
spatial and 8x higher mass resolution, we checked that 
the virial-infall times do change change significantly in 
these higher-resolution runs. This agrees with the res¬ 
olution tests in Garrison-Kimmel et al. (2014), which 
demonstrated the completeness of the satellite popula¬ 
tion at our mass range, M peak > 10 s Mg. 

A more significant concern is that ELVIS simulates 
only the gravitational dynamics of dark matter, and 
baryonic effects may change the survivability and stellar 
content of subhalos, in at least two ways. First, we as¬ 
sumed that all subhalos with M peak > 10 8 M 0 host lumi¬ 
nous galaxies, according to abundance matching against 
Af pea k, regardless of when the subhalos formed. We em¬ 
phasize that this approach is largely consistent with the 
observed mass function of satellites in the LG (Garrison- 
Kimmel et al. 2014), especially if one accounts for ob¬ 
servational incompleteness (Tollerud et al. 2008; Hargis 
et al. 2014). Furthermore, even if not all subhalos host lu¬ 
minous galaxies in reality, this would not bias our results 
if luminous galaxies populate subhalos largely stochasti¬ 
cally. However, some results from cosmological baryonic 
simulations (for example, Sawala et al. 2014) suggest that 
the subhalos that do host luminous galaxies are the ones 
that formed preferentially earlier, when they had deeper 
potential wells. If true, then this baryonic effect would 
shift the virial-infall times of surviving, luminous satel¬ 
lites to have occurred at higher redshifts. 

Second, the addition of the baryonic disk of the host 
MW/M31 galaxy can lead to more rapid disruption of 
satellites through tidal shocking or resonant stripping 
(Mayer et al. 2001; D’Onghia et al. 2010; Zolotov et al. 
2012). If a strong effect, then our dark-matter simula¬ 
tions would overestimate the survival lifetimes of satel¬ 
lites after infall, and the virial-infall times of surviving 
satellites would need to shift to lower redshifts. 

Thus, the combination of all potential baryonic effects 
could shift our virial-infall times in either direction, and 
future work should elucidate these trends with statistical 
samples of baryonic simulations. However, we do not 
expect that these baryonic effects would alter our results 
regarding reionization. 

After our submission, Koposov et al. (2015), Bechtol 
et al. (2015), and Martin et al. (2015) announced the 
discovery of multiple faint dwarf galaxies near the LMC. 


Our results strongly support the likelihood that many 
of these were satellites of the LMC prior to MW infall. 
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APPENDIX 

SATELLITES IN PAIRED VERSUS ISOLATED HOST 
HALOS 

While thus far we presented results using just the 
paired MW/M31 halos in ELVIS, here we compare our 
main results—virial-infall times and group preprocessed 
fractions—for satellites in the isolated versus paired 
MW/M31 halos. This comparison is useful for a num¬ 
ber of reasons. First, theoretically, we want to under¬ 
stand the degree to which the larger-scale environment 
around a MW/M31 halo influences the infall histories 
of its satellite population. Second, this comparison in¬ 
forms whether theoretical models need to consider sep¬ 
arately the satellite populations of paired versus iso¬ 
lated MW/M31 halos, for example, in order to under¬ 
stand satellites in the LG versus in isolated MW/M31- 
like galaxies in the local volume. 

Figure 8 shows ti r n s c t e infall (left) and 4 n ^ e /M3 i infaii ( ri § ht ) 
as a function of satellite M star or M peak , similar to Fig¬ 
ure 1, but for the isolated MW/M31 halos. For com¬ 
parison, dashed curves show the median values for the 
paired MW/M31 halos from Figure 1. While satellites in 
isolated MW/M31 halos first fell into a host halo slightly 
later, any such difference is small compared to the scat¬ 
ter. Thus, we conclude that the virial-infall times of 
satellites do not depend significantly on whether their 
current host halo is isolated or paired like in the LG. 

Similar to Figure 4, Figure 9 shows the fraction of all 
satellites at z = 0 that were a satellite in another host 
halo any time before falling into the MW/M31 halo (blue 
long dashed), or at the time of falling into the MW/M31 
halo (red short dashed), as a function of satellite M star or 
M peak ■ For comparison, the lighter dashed curves show 
the averages for the paired MW/M31 halos from Fig¬ 
ure 4 (left). Here, differences between isolated versus 
paired MW/M31 halos are stronger, such that low-mass 
satellites are more likely to have been preprocessed if 
they are in paired MW/M31 halos, at a level comparable 
to the halo-to-halo standard deviation (shaded region). 
This trend reverses slightly at higher mass, but here the 
difference is much smaller than the scatter, so we do not 
consider it significant. 

Most likely, the higher group-preprocessed fractions 
for satellites in the paired MW/M31 halos arises be¬ 
cause, as Garrison-Kimmel et al. (2014) noted, the paired 
halos have many more neighboring halos within a few 
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Figure 8. Same as Figure 1, but for satellite dwarf galaxies in the 24 isolated MW/M31 halos, which are matched in mass to the LG-like 
paired MW/M31 halos. For comparison, dashed curves show median values from the paired halos from Figure 1. We find no significant 
differences in the virial-infall times of satellites in the isolated versus paired halos. 
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Figure 9. Same as Figure 4, but for satellite dwarf galaxies in 
the isolated MW/M31 halos, which are matched in mass to the 
LG-like paired MW/M31 halos. For comparison, lighter dashed 
curves show averages from the paired halos from Figure 4. For 
low-mass satellites, those in the isolated halos are less likely to 
have been preprocessed in a group, at a level that is comparable to 
the halo-to-halo standard deviation. For high-mass satellites, any 
differences are well within the standard deviation. 

Mpc of them than the isolated halos, because the paired 
halos (almost by definition) reside in a preferentially 
higher-mass cosmic region. With more neighboring halos 
around, the satellites in the paired MW/M31 halos are 
more likely first to have fallen into a neighboring host 
halo. However, this difference in preprocessed fraction 
does not lead to a significant difference in infall times 
(Figure 9), so while group preprocessing is more preva¬ 
lent for satellites in the paired MW/M31 halos, the du¬ 
ration of this preprocessing is not longer. 

Finally, the most significant difference that we find be¬ 
tween the satellites in isolated versus paired MW/M31 
halos was in Figure 3: during the epoch of reionization 
(z > 6), the progenitors of the satellites in the isolated 
MW/M31 halos were much (~ 2x) closer to their nearest 
neighboring, more massive halo than those in the paired 


MW/M31 halos. This result may seem counterintuitive, 
given that the paired halos contain many more neighbor¬ 
ing halos at z = 0. However, we find that these structures 
were diluted over a much larger volume at 2 > 6 for the 
paired halos. Specifically, we randomly sub-sample all 
particles within i? v ir of each MW/M31 halo at z = 0 
and trace their locations back to z > 6, finding that the 
Lagrangian volume that contains all such particles was 
many (2 — 6) times larger for the paired MW/M31 ha¬ 
los, such that the satellite progenitors from the paired 
MW/M31 halos had fewer neighboring halos at a given 
distance at 2 > 6. 
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